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ABSTRACT 


A computer program was coded to obtain static aeroelastic effects 
on simple planforms through the use of subsonic lifting surface theory. 
The program was divided into two major technical areas, aerodynamic and 
structural, with matrix notation used to indicate the influence of each in 
the final aerodynamic loading distribution. Several typical stability 
derivatives were then obtained to make an accurate stability and control 


analysis on the desired planform. 
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I. INTRODUCTION 


In making a stability and control analysis on modern high perform- 
ance aircraft, the ability to make span load evaluations on the wing and 
tail with static aeroelastic effects taken into account is most important. 
It seemed desirable to develop a simple computer program utilizing 
Sument 'state-of-the-art' theories, and proper documentation, so that 
design groups might have readily usable methods at hand. 

Span load determinations have gone through a long line of 
evolution, right up to the lifting surface theories of today.  Prandtl's 
lifting line theory is considered as the earliest, which was suitable 
for simple planforms of high aspect ratio and unswept quarter chords 
operating in incompressible flow. Modifications of this method were 
made by Glauert, Multhopp, Weissinger, and others until the present 
lifting surface theories suitable for low (or high) aspect ratio plan- 
forms with swept quarter chords operating in compressible flow were 
attained. 

Two lifting surface concepts for subsonic calculations are: 

1. Define the surface pressure loading in the chordwise and span- 
wise directions by admissible sets of orthogonal functions with 
the magnitudes of the unknown coefficients being obtained by 
collocation techniques, Ref. 6. 

2 Define a grid of boxes over the wing, each of which contains 


a horseshoe vortex, Ref. 4. Then determine the strengths 





of the vortices by imposing the constraint of flow tangency at 

an equal number of control points. 

The approach chosen for this study was the second, which is 
usually denoted as a vortex lattice method and is thoroughly described 
by Belotserkovskii, Ref. 4. The major advantage to this approach 
is the convenient use of matrix methods for taking the static aero- 
elastic effects into account. 

The method of Gray and | Ref. 1, used initially on the 
Boeing B-52 aircraft, was carefully reviewed as background for this 
study. Their analysis, which was restricted to aspect ratio wings of 
greater than five, employed the Weissinger L-Method to account for 
subsonic aerodynamics. The advantage of the lifting surface theory 
used for this study is that wings of lower aspect ratio may be ana- 
lyzed and the ability to more accurately predict chordwise (as well 
as spanwise) distribution of aerodynamic loading is provided. The 
method of Ref. 1 did not allow for aerodynamic loadings other than 
additional бе! whereas the method of this study is not restricted 
to loadings with the aerodynamic centers at the local quarter chord. 

To account for compressibility in subsonic flow the Prandtl- 


Glauert rule was used, Ref. 2, sometimes referred to as planform 


Daa tere type pressure loadings are defined as the air loadings 
due to angle of attack only. Generally they have a cotangent type chord- 
wise pressure distribution in subsonic flow with a logarithmic type singu- 
larity at the leading edge, and vanish at the trailing edge to satisfy the 
Kutta condition. 
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distortion. The aeroelastic effects will depend upon the original wing 
structure and will vary proportionally with the free stream dynamic 
pressure. The span load evaluation at zero dynamic pressure corres- 
ponds to the rigid wing results. Symmetric span loads will be appli- 
cable to longitudinal stability and control analyses, with anti-symmetric 
to lateral-directional. 

The program has been restricted to simple planforms with 
straight taper, and the influence of fuselage interference was not 
taken into account. This was not an oversight but deliberate since 
it was felt appropriate at this time to obtain a straight-forward docu- 
mentation on a simple wing situation, with applications to more com- 
plex planforms and interference effects as a logical extension of this 
study. 

Ihe approach used in the structural area was restricted to a 
basic beam-like wing structure, corresponding to modelling the wing 
by an elastic axis with variable stiffness properties along the span 
for both bending and torsional moments. An influence coefficient con- 
cept for plate-type lifting surface might be a more general approach, 
but this was not felt warranted without a specific application in mind. 

The purpose of this study was to develop and document a simple 
approach for static aeroelastic analysis and airframe stability and 
control considerations. This documentation is a necessary prelude to 
more serious considerations on actual configurations such as new 


current Navy aircraft like the F-14 and S-3. 
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П. METHODS USED 


A. GENERAL 

In a static aeroelastic analysis there are two distinct technical 
areas requiring attention. Both the aerodynamics of the problem and 
the structural aspect must be considered. 

With the free-stream subsonic Mach number and dynamic pressure 


given, the aerodynamic loadings of the rigid wing are given as: 


[a] {St} = fox} (.01) 


where [ A] 


Aerodynamic influence coefficient matrix at the given 
Mach number. 


A = incremental loadings over the wing. 


E 


|} 


rigid wing input angle of attack, specified at 
the control stations. 


Structural deformations on the wing due to the aerodynamic load- 


ings result in a change in the angle of attack. 


q isi {Sf} = [&s} (.02) 
where q 


[s] 


(0%) = the angular deflections at the specified control 
points on the wing due to aerodynamic loadings. 


the free stream dynamic pressure, scalar constant. 


the structural influence coefficient matrix. 
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Now the complete solution to the aeroelastic problem is given by: 


йа1- «tsi {H} = {x — со 


al- [TAT- tsi] [e д) MEE 


B. AERODYNAMICS OF THE PROBLEM 

The basic approach used in this development is similar to the 
lifting surface analysis as described by Belotserkovskii Ref. 4. 
With the concept of a vortex lattice array of horseshoe vortices the 


wing is covered with N boxes where: 


N = NCxNS------------ (.04) 
and NC = Number of chordwise station 
NS = Number of spanwise station 


In a manner similar to lifting line theory, as developed by Prandtl, 
Glauert, Multhopp, and Weissinger (Refs. 2and 3), the swept bound 
horseshoe vortices were placed at the local quarter chord of each box 
and the flow tangency requirement imposed at the local three-quarter 
chord control point of each box. Figure l1 shows the geometry of the 
problem, with A and B denoting the left- and right-hand corners of 
the horsehoe vortex element and C the control point. 

As a consequence of the Helmholtz vorticity theorems (Ref. 2), 
and because inviscid flow was assumed, each horseshoe vortex ele- 


ment has constant strength, once established. The Biot-Savart Law 
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Figure 1. Horseshoe vortex lattice geometry. 





Figure 2. Finite vortex element, between A and B, 
of strength 
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(Ref. 2) then defines the velocity induced at a control station C (see 


fig. 2), for a vortex filament of finite length as: 


[5 
мш = C > coe) er: (.05a) 


and the downwash angle as: 


Aw _ Г ) 
E iG - ) (cos x + cos ë Pm (.05b) 


If the vortex element is expanded to plus and minus infinity, the 
well known expression for downwash angle in two-dimensional flow 
Aw 
results: = Fe SN ‚06 
p" NC To V al (-06) 
Eiuee  cosc = cose = 1.0 
As previously stated, the flow tangency requirement was applied 


at the three-quarter chord control point, and the flow was required to 


be parallel to the zero lift line at that point. Hence, 


О (.07) 
E = f 


Justification for the use of a local three-quarter chord control 
point in each vortex lattice box has not been rigorously established. 
The concept of the "rearward" aero-center at the wing three-quarter 
chord control point, from steady and unsteady two-dimensional wing 
aerodynamics, appears to be the genesis of this usage. At present 
the modelling is considered satisfactory because the theoretical load- 
ings are in accord with experimental results. It is mot clear as of 
this writing why the preservation of flow tangency at local vortex- 


box three-quarter-chord control points preserves the Kutta condition 
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at the wing trailing edge, but the calculated values of local vortex 
box loadings do show the preservation of this trait. (Ref. Appendix D) 
LO ; 
The value of the downwash angle a varies at each control 
E 
station on the wing, and there are contributions by the bound and 
trailing vortices from all of the wing stations. Hence, the flow 
tangency requirement may be expressed as a column vector as follows: 
wy = + С} mE 52. dos) 
{ Я И ГАУ q 


4 


* 
where the matrix [А] is the aerodynamic influence-coefficient matrix 
as developed in Appendix B. 
For the vortex lattice representation with N horseshoe vortices, 


the flow tangency requirement was given by: 


(y= (aI = aq 


=— 


er (2) = la] (e u... (.09) 


There was a temptation to view the vortex strength at each 
control box as the pressure differential of that box. This, however, 
is not the case. Again as a result of the Helmholtz vorticity theorems 
(Ref. 2), the vortex strengths of the individual boxes may be summed 
in the chordwise direction to obtain the total vorticity at each span- 
wise station. Hence the vortex strength of a section may be "lumped" 


into one vortex equal to the sum of the individual vortices. (See Fig. 3.) 


* The development of the Iu matrix for spanwise stations 
only may be seen in Ref. 1l. 
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Figure 3. Contour encircling all ls located at the 
individual boxes. 


a) Field X, Y, Z b)" Field'Xo, Yra 


Figure 4. The effects of planform distortion. 
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From the total vortex strength at a chordal station, the section lift 
NC 
loadings can be obtained, i.e., a= > АГ. с 10) 


t= A 
and by application of the Kutta-Joukowski law, 


(5) = 2 а er) 


The section moment may be obtained by taking the moments about a 
reference station of the individual bound vortex elements, inasmuch 
as the bound vortex elements represent the local air loading of each 
box by the Kutta-Joukowski law. 

A two-dimensional study was made to compare the "pressure 
distribution" imposed by the vortices and the dem type pressure 
distribution as prescribed in theory (Ref. 2). 

Assuming that the chordwise vortex elements described a chord- 
wise pressure distribution, it appeared feasible to use a chordwise 
collocation technique for obtaining section lift coefficients and moment 


coefficients, C, € men. given by 


C," er A, a m,, € c S uS 


mea 
A E. AS р 
Since x, > c7 =! cos ©) 
апа eX <= AA сот © BURN Sin nO 
rae o xn 
Y e nel 


it appeared that a solution for the Fourier coefficients could be found 


in the following relationship: 
cov @, SING, SINZO, --- SINN-DO, 
| | d TER 
| Le) = | | | IT ben [А] 
q | | үм Von I 
| | | | 
сот e S\N On SIN 26, n SINAD on 
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However, the existence of the A,'s assumes knowledge of the pressure 
distribution over the wing, and correct results were not a, 
since the vortex strengths did not represent pressure loadings. 

A check was made using two-dimensional theory to bear out 
these conclusions, and it was found that even for ten chordwise sta- 
tions, the pressure distribution was incorrect, and the aerodynamic 
center was not located at the one-quarter chord. The leading edge 
singularity was not properly represented and the first collocation 
point differed from theory significantly. This check can be seen in 
Appendix D. 

In light of these results, a straight-forward summation was used 


to obtain the section lift coefficients , i.e., 


NC A 
Cc ZI - — — = (12) 
c= 


©. COMPRESSIBILITY CORRECTIONS FOR SUBSONIC FLOW 

The compressibility corrections were obtained by the Prandtl- 
Glauert planform distortion method (Ref. 2). The applicable flow field 
equation for the velocity potential was given by 

2 

(1- M2) ш + B dE = о MENSES 
where the subscript represents second partial derivative. By use of 
the Prandtl-Glauert transformation, X = XN I- MP = = Жш 


this reduced equation (.13) to: 


$ + Фф +o =O "m (15) 
o gy 22 


Xo X 
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The effect of the Prandtl-Glauert transformation was to stretch all of 


| 
the x-coordinates by the factor NUT . The effects can be 


2 
l- MS 
seen in Fig. 4. 
The net effects of the planform distortion were as follows: 
1) the aspect ratio was decreased; 2) the leading edge sweep was 


increased; 3) the taper ratio (Смс, ) was unchanged; and 4) most 


importantly, the downwash distribution {2} was unchanged, i.e., 


Q- 


e e. 


| ea” 


Since the downwash angle distribution, [2| , is unchanged in 


Dm ae 


us o is M (.16) 


going from the distorted wing to the original wing, the aerodynamic 
bsec: AU š ; DS 
loading distribution, q! , will be the same on both wings. This is 
an important concept in that the aerodynamic loading obtained from 
the distorted planform under incompressible flow will be the same as 
the aerodynamic loading on the original planform under compressible 
flow. Hence the aerodynamic loading is an invariant when going from 
the compressible undistorted planform to the incompressible distorted 


planform. Б 


р. STRUCTURAL ASPECTS OF THE PROBLEM 

Because wings are flexible or elastic, there are always changes 
in angle of attack due to the deflections of the structure. There are 
two deflection modes present, torsion and bending, both affecting the 
changes in angle of attack distribution on the wing. The effects of 


these two modes were called static aero-elastic effects and have been 
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incorporated into the program by the simple elastic-axis beam 
representation. 
The final angle of attack distribution [o] was made up of 


two components CA , and EA , 1.€. 


СЯ = CA i { cx, } Оо---- (.17) 


Now the (o) as shown in Ref. 1, has been determined to be: 


{as} = sy 2 Em (.02) 


See Appendix C for the development of the structural influence 


coefficient matrix [s] : 


Е. TOTAL SOLUTION TO THE PROBLEM 
By combining equations .12, .17, and .02, the following rela- 


tionship was formed: 


[a = als) + leq] 
E [La] - q[s]] Lo = [о | A (.18) 


The closed form solution for aerodynamic loading distribution 


may be given as: 


a “sp m te] " 3 (.03b) 


This was the final equation used in the program to obtain a chordwise 


and a spanwise loading distribution for a given surface. 
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F. CALCULATION OF STABILITY DERIVATIVES 
Five stability derivatives were investigated to check the validity 


of the program. These were: 


Co. S CS M © 


© e Ce Sa 


The total wing Сү, was obtained by a spanwise collocation 


technique, using the section lift coefficients. 


+1 
L 
25-6 eui) 4) 4% = zn 

where ne; dy=-sine de 

E Y = ei B_sın n9 
SN. 

Th eo . 

en = 2; B, ew no sino de 


T 
ЕК 
B, j i nO sin O дө = 7 9 


Since by orthogonality 


Tr TY coe nel 
S sw мо мо ао = } 2 
O a nyi 
Therefore 
TX B, | 
= eee e as з "20 
C. A Cane e 


Since the assumption of linearized wing theory was made, 


E = e due [O fo. input of 1.0 rad. uniformly over the 
E NE 
t. 


whole wing, i.e. C ^а FF 
[x= to] 
The Cm, 5 were obtained by a straight-forward integration 


technique. In effect, the section pitching moment at the quarter chord 
NC 
om 2 Ar 
was found to be: ( QA = Cm mn = p (RK) 
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Hence, total wing pitching moment is: 


C uo E UB olx) 


Again C was evaluated by: Cm = Ne 


„= аш (.21) 


where AC. is pitching moment change about the wing a point 


due to {a} = 1.0 rad. input. 


With ÉL and C m calculated, the Cm. was easily found 


by the following relationship: 





dc C 
С = м = =m —— š 


From the anti-symmetric solution of the Aero Matrix (see 
Appendix B), certain rolling moment derivatives were extracted. The 
rolling moment coefficients were obtained by a collocation technique, 


similar to that used to obtain the total wing er 
+I 


= = -x YQ) nde | 
" = E Ще "nO coso sinode 


=\ 


LS е ° u20do 2 -L-B 
" = Tm 2 B., sin M ©\ lo e 


TY 
TV, For n=2 
өдө = 2 
where 1 Sin NO SINZ | ТЕ 


° 


ER A: B2 _ (.23) 
Therefore C= Ve e — = = 
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The roll damping derivative, Ce, , was then obtained by placing 


the anti-symmetric angle of attack input on the wing, as shown in 


pb 
Figure 5. At the tip, the angle of attack - Z V = ] rad., and 


Co was given by: 


є „3% AG 
Ar ME CASA 


P 


A (.24) 


The aileron control effectiveness, em , was found by defining 

the location of two of the boxes to coincide with the location of a 

typical set of ailerons, see Fig. 6. By inputing zero angle of attack 

everywhere on the wing except at the aileron locations and giving 

them 1.0 rad. angle of attack each, positive on one side of the wing 

and negative on the other, a rolling moment was produced. e Was 
а. 

l given by: 

EN, 


(= = EL. l2. o ANS) 
Ay So {| OL, DUE To Sa} 
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Figure 5. Angle of attack distribution for a constant 
rolling velocity P. 





Figure 6. Ailerons as defined by horseshoe vortices. 
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ПІ. RESULTS AND DISCUSSION 


A. LONGITUDINAL 
The results from the symmetric portion of the program were 


checked by using longitudinal stability derivatives (i.e. lift-curve 





ОСЕ 
а СЕ 


that changes in leading edge sweep and Mach number are handled 


slope er , moment-curve slope Cm, , and To insure 
correctly in the program, checks were made on C. and eu. 
These checks were made for a rigid wing with an aspect ratio of 6.0 
and a taper ratio of 0.3. 

Figure 7 shows the effects of leading edge sweep on lift-curve 
slope C. , and o cures slope Emo . These curves were 
obtained by inputing a "do loop" in the program and varying the lead- 


ing edge sweep ji from 0 to 40 degrees. These results agreed 


( GM =o 
VEME cos Aey 


quite well with theory, i.e. ( l \ = sag (26) 
AA 
where the Mach number is held constant and the sweep varied. 

Figure 8 shows the effects of the compressibility corrections in 
the program on C. and Cm . Again a "do loop" was used, 
stepping through the program for Mach numbers of 0.0 to 0.9. These 
results also agreed with theory, i.e. equation (.26), where the sweep 


is held constant and the Mach number varied. 
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Figure 7.  Lift-curve slope plotted as a function of leading edge 
sweep at a Mach number of 0.5. 
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Figure 8.  Lift-curve slope and moment-curve slope plotted as 
a function of Mach number with a leading edge sweep 
or 20.0 deg. 
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It should be noted that the Emo variation, with increasing 


leading edge sweep, is in a direction to shift the aerodynamic center 
ee 
dC. 


about the wing C/A corresponds to the distance, in wing m.a.c. 


aft. Figure 9 shows this more clearly since the value of 





lengths, that the aero center is from the reference axis. A positive 


dC, 
ver 


C/a reference point. 





value of indicates that the aero center is forward of the 


B. LATERAL 

In order to document the anti-symmetric portion of the program, 
two lateral stability derivatives were investigated, roll damping Sn. 
and roll due to aileron deflection Ce. . The aeroelastic effects were 

a 
checked at the same time by choosing an example from Bisplinghoff, 
Ashley, and Halfman, Ref. 3. A typical jet transport with a straight 
taper and straight elastic axis defined at 3576 c is E. The ÄR for the 
wing was 6.15 and the taper ratio was 0.444, The bending and torsional 
stiffness properties of the wing were taken from Ref. 3, Fig. 2-29; then 
they were scaled to our model by (EI) = 2 (ET), and reproduced 
here in Fig. 12. 

The checks were made at a constant Mach number while varying 
the dynamic pressure (q in p.s.i.) from 0.0 (rigid wing) to 18.0, corres- 
ponding to varying the altitude. The effect of aeroelasticity upon roll 
damping for this wing, with an unswept quarter chord, is to increase the 
damping with an increase-inq. Note that the sign on roll damping is 


negative since the moment that is developed is opposite to the direction 
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Figure 9. 


AR=6.0 
TR=0.3 
M=0.5 


10 20 30 40 
IA T ADISIGU 
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Static margin indicator plotted as a function of leading 
edge sweep at a Mach number of 0.5. 
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ofroll. Figure ll shows the roll damping derivative plotted as a function 
of q. The results of aeroelasticity upon the rolling moment due to aileron 
deflection can be seen in Fig. 10. Note that the change in sign corres- 
ponds to the control reversal dynamic pressure, Чг. The q, Was found 
to be 11.8 p.s.i., while it was shown to be 11.39 p.s.i. in Ref. 3, by 


another method; this is a difference of approximately 4%. 
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Figure 10. Rolling moment due to aileron deflection plotted as 
a function of free stream dynamic pressure at a 
Mach number of 0.5, 
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Figure 11. Damping in roll derivative plotted as a function of 
free stream dynamic pressure at a Mach number 
ОО. 5. 
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Figure 12.  Bending and torsion stiffness curves. 





IV. CONCLUSIONS AND RECOMMENDATIONS 


This computer program can provide a 'limited' capability to make 
stability and control analyses of modern Navy aircraft. The program, as 
it now stands, is limited both to wings structurally compatible to the 
elastic-axis concept and to simple planforms. However, slight modifi- 
cations could easily handle arbitrarily shaped wings and even a variable 
sweep configuration, such as the F-14. 

The study is not complete at this point; a logical addition to the 
program would be to add the influence of the body radius to the Aero sub- 
routine. This would improve the value of the program. If felt warranted, 
plate theory could be used for part or all of the structures matrix. 


A follow-on study now seems appropriate to make the program more 


attractive and useful to design groups. 
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APPENDIX A. ASSUMPTIONS 


In writing this program, many theories with basic underlying 


assumptions were used. Most of them are associated with thin airfoil 


theory, or the structure of the wing. These assumptions, most of them 


listed by Gray and Schenk in Ref. 1, are as follows: 


1: 


Potential flow about the wing with no boundary layer effects , 
i.e. inviscid, and without compression shocks, i.e. below 
transonic mach numbers. 

The thickness ratio of the wing is small. 

The lift-curve-slope, @ = CTX, i.e., two-dimensional 
Changes in camber due to bending and torsion of the wing 
are neglected. 

The changes in angle of attack due to structural deformations 
are uniform along each spanwise station. 

The elastic deformation of the control surfaces i the same 
as that of the adjoining wing. 


The angles of structural deformation are small such that 


lle 


SIN © = tan OF © cos GC =1.0 


> 


Drag-load effects to deformation of the wing were neglected. 
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APPENDIX B. DERIVATION OF THE AERO SUBROUTINE 


The basic development of the downwash matrix or aerodynamic 
influence matrix was done by Gray and Schenk in Ref. 1. An expansion 
was made here to include swept-bound vortices and a chordwise, as well 
as spanwise, distribution of horseshoe vortices. 

The lift or aero loading distribution over an airfoil can be charac- 
terized by a system of horseshoe vortices, each of these vortices con- 
tributing to the upwash or downwash at points all over the wing. Figure 
Bl shows two horseshoe vortices combined to make a wing. Conventional 
right hand rule indicates the sense of circulation associated with each 
horseshoe vortex. Quite obviously, the more horseshoe vortices used to 
approximate a wing, the more accurate the lift distribution obtained. 

With the wing divided up into NC X NS stations, the horseshoe 
vortices were placed at each quarter chord of the box, and the horseshoe 
vortex lattice, as shown in Fig. 1, was arrived at. Now а boundary 
conditions or flow tangency requirements, for simplicity of calculations А 
were imposed at the three-quarter chord Í for each station or box. The 
boundary condition implies that the downwash angle ( =) at each cons 
trol station is equal to the geometric angle of attack. The downwash vel- 
Ocity UJ at each station is due to the sum of the downwash velocities in- 


duced by each horseshoe vortex over the entire wing. 


ЫЕ chord for each box is equal to the local wing chord divided 
by NC. 
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Figure Bl. Simplified horseshoe vortices system for 
wing lift distribution. 
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Figure B2. Geometry for downwash velocity calculation. 
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Given one horseshoe vortex and one control point on the same 


wing, as shown in Fig. B2, the downwash velocity Aus was calculated 


as follows: 

Учса Rey y OS 
= (ХИ SING = (X.-Xe)/Re 
cos @= (у-Ч„)/ Ед cos) = (u-- Ya)/ Re 


0 ъ= \ (О + (Xe- XA)” 


SIN Y= (Xe EN cos Y= Cys- Ya) / 9 дв 
х= е-х € B=0-Y% 
SIN X = чм Ө cos X — cos Ə SIN Y 


COS CX — cos@ cosxX¥ + SNO sN Y 


cos B= cosd cosy + sin ф ain У 


The left hand vortex contribution to the downwash velocity was given by: 


T (l+ sim o) 


AW, = БЕ (В.01) 
Ar ( Yor Ya) 
The right hand vortex contribution was given by: 
NS 2-Е (1 емф) 
p ML uice > (В.02) 
Дт ( Yo.” us 
The bound vortex contribution was given by: 
(^ (cosa- cosB) 
NuUa m 32 == (B. 03) 


An Ry sinc FT 
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If the control point happened to be on a straight line extension of 
the bound vortex, both numerator and denominator would go to zero 
simultaneously. The contribution from the bound vortex would be zero 
and a check was placed into the program and if this situation were to 
develop, MAW, would be set to zero. 

The total contribution to the downwash velocity due to the one 


horseshoe vortex on the same wing was: 
AW=AWS + Аш + A tos 


Or ^ uJ = dui (14 SIN Q) " ( |+ SIN d) — (B.04) 
Arv * (y. - 4) (Y-Ye) Rasa 


The downwash velocity contribution at control points all over the 
wing, due to a horseshoe vortex system, was calculated in a similar 


manner, These can be thought of in matrix format as follows: 


W, AUD + AW. + d oe Ао м 

Wwa NO, t AW . . . + AU 
= l 

| | 

Wu AUD + А04 + Š s ye Xn + iU 


where the subscripts indicate the downwash at station i due to a horse- 
shoe vortex at station j, etc. Now by factoring a n from each term 


in the proper manner: 
{= АГ} 


then the following relationship follows: 
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DE [А1 2) | wu | (B.05) 
or (3) = LA] {St} M (B.06) 
Now {= [о] "ы Гос = ТАТА BO) 


The A matrix is the aero-influence matrix as used in the program. 


The computations are carried out for contributions from each wing, i.e 
ue NI 


For symmetric loads: L AS\= | Agl *| Де | (B. 03a, b) 
For anti-symmetric loads: т А N) = L A Е А .| 
Then LAS Las OR ANM 


Also contained in the Aero subroutine is the Prandtl-Glauert plan- 


RH. SIDE 


form distortion correction for subsonic compressibility. In essence, all 
T E 
of the x distances are stretched by a factor of (1- м5) - , defined as 
e Е Therefore throughout the subroutine 5(Bx,y) replaces f(xy). 
e. Qc x4 instead of (X,- X4)  ). 

The 'lifting-surface' representation of a wing has been written up 
in FORTRAN IV for use with the IBM 360-67 computer. The program was 
designed to hand]e planforms of low to high aspect ratio, swept leading 
edges, and a linear taper ratio. The choice of the number of stations on 


the wing, chordwise and spanwise, was left arbitrary for easy compati- 


bility with most wings and for almost any flap and aileron arrangement. 


4] 





So as to make the program more flexible, subroutines were used 
for each phase of the calculations. Subroutine Geom is used for locating 
the horseshoe vortices and their control points; subroutine Aero is used 
for К ана the aerodynamic influence (or downwash) matrix; and sub- 
routine Struct is used for obtaining the structural influence coefficient 
matrix. With the subroutine packages set up as they were, changes to 
the program can be easily handled simply by making small changes in the 
input or the main program. 

Sample calculations were made for arbitrary wing geometries. The 


results of these calculations can be seen in the main text. 


The desired results were obtained in these calculations. The 
corrections for subsonic compressibility, sweepback, and static aero- 


elastic effects have been checked and do give the desired accuracy. 
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APPENDIX C. DERIVATION OF STRUCTURES SUBROUTINE 


As in Appendix B, the [s] , Structures matrix or elasticity matrix, 
was developed by Gray and Schenk in Ref. 1. An expansion to fit 
‘lifting-surface' theory vice 'lifting-line' theory was made for this 
program. 

The problem was basically stated as: Given a loading over a wing, 
find the deformation imparted to the wing and apparent change in angle 


of attack. From the previous subroutine, 


м] 
o = Lalis ONUS. (B.07) 
The following relationship was obtained from this subroutine: 


fa} = m n (.02) 


The assumption has been made that the wing acts like a beam with 
the elastic axis representation. (See Refs. land 4.) The beam bending 
and torsional moments have been resolved about the elastic axis at each 
spanwise station of the wing. Although the elastic axis may be arbitrary, 
it was thought of as a straight line for simplification here (see Fig. Cl). 

The wing is divided evenly in the spanwise direction with each of 
the chordal strips having width 2h. Therefore for the wing with an 
assumed unit value of semispan, h = 1/2 (NS) EA den (C.01) 


Positive direction for torsion and bending moments are shown in Fig. Cl 


relative to the elastic axis orientation. 
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Figure Cl. Elastic axis representation of the wing with 
the torques and moments at the station elastic 
axis points. 
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The span loads per unit span length are defined as М for the 
N vortex lattice boxes in a numbering sequence starting from N = lat 
the most forward, inboard box proceeding aft to the trailing edge at 
N = NC, then stepping to the next outboard set of boxes at the leading . 
edge box for N= NC+1. Continue in this manner until the Nth box at 
the trailing edge of the tip is reached. 

The individual lift force contribution towards the bending of the kth 


wing station may be expressed as: 
‘Смс (Uy 


Dig = (45) (2- 3,4) qh A o PT (C.02) 
anos  ! 


The presence of the Kroecker delta takes into account that only one-half 
of the air loading at the kth station contributes to the bending (and torsion) 
moment at the kth station. 

Therefore the bending moment at the kth station may be expressed 


as; 


NS 
o | 
M, = AL AZ AL, EWR) ___ (0.0) 


or in matrix notation as: 
(m] = q Lam1(G] EL s 0M) 
where {m} = N5xt column vector 
l. = N« 1! column vector 


CAM) = NS«N rectangular matrix, and 
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cs ———— € — 0 | —— ee 


[Ам] = © 


The AM matrix was partitioned into NS x NS elements where each ele- 


ment is a (1 x NC) row matrix made up of coristants as indicated, i.e. 


ВЕ | 


a = lau... Vi 
Bc = 


The torsion moment at the kth station may be expressed as: 
= A | 
Tee aln (.05) 


| | 
Ku: Mine AX PAX 


where 


— 
 — 


= n = 
ти ш иг zen 


and 


== m E ( Key); 


where 
I$ 4NS. € 15 < Мамсхмѕ 


The zero partition elements to the left of the lead diagonal correspond 


to inboard loadings that do not contribute to the moments at stations 


outboard. 
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The torsional and bending moments as developed here were perpen- 
dicular and parallel respectively to the streamwise flow at the elastic 
axis reference points. In order to reorient them relative to the elastic 


axis, they were put through the following rotation matrix: 


M' Е соѕ5 Лед -SiN Aea M 


— Е CANO) 
T ОСЕ s se s 


Now the moment and torque distributions have been resolved to be the 


following M cuu 


[w] = а [es А Tem AcadaT tt 
(y) Е | їз М AM |+ [cos Aes ara) 


The change in angle of attack in the streamwise direction was ob- 


-- (C.07) 


tained by the relationship (see Ref. 1): 


n n : 
- \mMds ,\ tTds (C .08) 
= ERI | GJ 


where m equals the beam bending moment per unit pitching moment at 


the station n , and t equals the torsion about the elastic axis per 


2% 


== Е CG . 09 
со» Лед ( 


unit pitching moment. ds > 


The above integral equation was changed to an algebraic equa- 


tion by summing the contributions at each station. 


= Quis + 2 mu Meta er E ) 


sy (ED, cost, (ЕЛ), соѕЛ, 
иво T f. A Re ) (С.10) 


(G3),cos4, (GJ), cos A, 
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The influence coefficients mij and tij were easily defined by the geom- 


etry of the problem and found to be as follows: 


-sin A; where the point i is outboard of the point j 
m | 
0 where the point i is at or inboard of the point j 
cosA; for station i outboard of station j 


and tij = 
0 for station i at or inboard of station j 


In matrix form the mis and tu's were written as follows: 


= ux des = 


2 mó - SIN A, "Ыы МС сс 
Ó - sin Ar -SiN Aa Cee eg 
Е j A 
6 О ~SinAs _ __ 
! | E 
| ! : 
cos, cos, ca. = 
2 
BE O со5Ла coshz ____ 
2 
a со Лз 
| 
! 


Now 


[o х. = [mf — “ м} 1107 (С.11) 
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Thus the desired form was achieved and the S matrix was: 


le Lm ][eos/ [AM] = ew A] ol 
m mu a a] fan A JAM] +feosA] А | 3220 


The final resultant matrix equation was: 
= Al | 
{о} = 4 15109 о 


The computer program is normalized such that the wing semispan 
has unit length. However, an actual wing will have a finite wing span 
of say b/2 inches, will operate at a dynamic pressure of q p.s.i., and 
will have stiffness properties of EI and GJ in units of eine е 

Recognize that in a scaling analysis the structural twist of the 
model will equal the structural twist of the prototype and that they will 
operate at the same dynamic pressure. However, the aerodynamic load- 


ings will differ by a factor of А 


(94). b | 
where À = ea ЕЕ = ( a, since (2) == | (C.13) 


m 


(9 - 1, ca 


Therefore by equation (.02) we can scale 


LS ae Sle 
but analysis of equation (C.12) reveals that {sloc WEL , 9/67 
[S]. С [Ee]. 
so Nu. Fir]. 


UU $e 2 Ges X (ED | 


[ 
о --- (C.15) 
"e ПЕ = x (GJ) | 
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APPENDIX D. TWO-DIMENSIONAL WING PRESSURE 
LOADING CONVERGENCE CHECK 


In this appendix, a check was made to ascertain the validity of 
viewing the vortex strengths of the control boxes as pressure differentials. 
A comparison was made for the Cp S obtained by this assumption, using 
a collocation technique for obtaining a chordwise pressure distribution, 
and the theoretical Cp $ with a cotangent SETS pressure distribution. 

A two-dimensional section with unit chord was chosen for this 
check. The section was divided into N equal chordwise stations with 
the vortices placed at the quarter-chord of each station and the control 
points at the three-quarter chord of each station (see Fig. D1). 


p 4 


Since (8 ® represent the small vortex strengths concentrated 


) ) 


at the quarter chord for each station, the total vortex strength was given 





by: 
Al; EVR 
N Al,=eVT: AP LLLI 
И DEN e £ £ Ne C/N 
P; 3 
and AC, = EU Um (E mem. (D.01) 
t 1⁄2 uS Q/w) ^ V 
The principle of an 5 vortex (Ref. 2) was used for the down- 
(ç) 
wash angles, i.e. y = e or 
о 
(l= = Le liz т дЕ (р. 02) 
where v ОС at the three-quarter chord control point, and 
e (D.03) 
C 
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Figure D1. Geometry for two-dimensional study of 
chordwise pressure distribution. 
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with ХА; = = j- => , the local quarter chord point, 
and 55 + bo = 2.25 , the local three-quarter chord point. 
t 


Hee — (AC,|] 2 Anu lel (e) ——_- 0.04 


with K = 1.0 radian for each station. 


Now the theoretical A G for two dimensions was given by 


(Ref. 2): ANC S = 4 (1 + соф Ө} 





— NUES . (D.05) 
using the coordinate transformation: 

š = 5 ( 1- cos ©) 

E o 2(%] пиу 


Since the check was made at & = 1.0 radian and a unit chord, the 





Сус obtained will be equivalent to Cry and should be approximately 
Cm. e 
2M The aerodynamic center given by: AC = Ge (D.07) 
1 


should be at 0.25c or just 0.25 inthis case. 


The tabulation in Fig. D2 clearly shows that the correct results 


were not obtained by the method as described herein. 
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SUBROUTINE ELU TRI-DIAGONALIZES THE INPUT SQUARE MATRIX FOR SOLUTION. 
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